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Abstract 


Recent field and analytical studies of what has been 
traditionally called “Paliza Canyon obsidian” in the 
archaeological vernacular show it to be Bearhead 
Rhyolite that is part of the Late Tertiary (Neogene) 
Keres Group of the Jemez Mountains, northern New 
Mexico. The geological origin of all other archae- 
ological obsidian sources in the Jemez Mountains 
have been reported and are well documented in 
the literature. But the so-called “Paliza Canyon” 
source, important as a toolstone to Pueblo Revolt 
Colonial period occupants of the Jemez Mountains 
area and present in regional archaeological contexts 
throughout prehistory, had remained unlocated and 
undocumented. The Bearhead Rhyolite origin for the 
“Paliza Canyon” obsidian (which we suggest should 
now be named “Bearhead Rhyolite”) solves this 
ambiguity and provides more precise geological and 
geographical data for archaeological obsidian source 
provenance in the region. 


Introduction 


An important type of obsidian used by Native Americans, 
traditionally called the “Paliza Canyon” obsidian 
(Nelson 1984; Baugh & Nelson 1987; Wolfman 1994; 
Glascock et al., 1999; Shackley 2005), is widespread in 
the archaeological record in New Mexico. Church (2000) 
and Shackley (2012a) have both reported the presence 
of Paliza Canyon obsidian in Rio Grande Quaternary 
alluvium as far south as Las Cruces, New Mexico. It also 
occurs commonly in Pueblo Revolt period archaeological 
sites (late 17th century) in the southern Jemez Mountains 
of northern New Mexico, including near Jemez Springs 
(Fig. 1; Liebmann 2012; Shackley 2009a, 2012b). Using 
x-ray fluorescence spectrometric (XRF) analyses, these 
archaeological samples have been chemically correlated 
to a “Paliza Canyon” geologic source defined in the Baugh 
and Nelson (1987) study. XRF is the preferred method for 
chemically characterizing archaeological obsidian because 
it is relatively inexpensive and can be accomplished non- 
destructively (Shackley 2005, 2011). Determining the 
eruptive location of this unique archaeological obsidian 
is crucial for understanding exchange, migration, and 
social networks during the Pueblo Revolt period (A.D. 
1680-1692), as well as all prehistoric through Colonial 
periods in New Mexico (Liebmann 2012; Shackley 2005). 

The original four samples of “Paliza Canyon” obsidian 
were collected in a portion of Paliza Canyon, a drainage 
located in the southern Jemez Mountains (Figs. 1-2, near 
our sample locality 100915-2; Nelson 1984; Baugh and 
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Nelson 1987). It is important to note that these samples 
were from late Quaternary, valley bottom alluvium and 
thus geologically reworked. The goal of our study was to 
locate the primary geologic source of these samples. 

The original “Paliza Canyon” source was characterized 
as a small-nodule Neogene source, based on the age 
of adjacent volcanic rocks and small nodule sizes (as 
opposed to large nodule Quaternary sources, such as 
Valles Rhyolite obsidian). The small nodule sizes likely 
reflect hydration and perlitization of volcanic glass over 
relatively long periods of geologic time. Known Neogene 
and Quaternary sources of archaeological obsidian in the 
Jemez Mountains have been well dated and chemically, 
stratigraphically, and isotopically studied for a number 
of years. The geologic context of the “Paliza Canyon” 
obsidian has been more difficult to decipher because of 
complexities in its geochemistry and geologic setting 
(compounded by its relatively old age), especially when 
compared with most other Jemez Mountains sources 
(Gardner et al., 1986; Self et al., 1986; Spell and Harrison 
1993; Ellisor et al., 1996; Goff et al., 1996; Kues et al., 
2007; Kelley et al., 2013). 

Relatively recent geologic mapping and sampling 
efforts in the southern Jemez Mountains, combined 
with XRF analyses of several samples, indicate that 
the primary source for the “Paliza Canyon obsidian” 
is only a few kilometers north of the original proposed 
source area. After providing a geologic background, we 
use geochemical data to compare the “Paliza Canyon 
obsidian” with possible primary geologic sources in the 
area. We argue that the best candidate for a primary 
source is a geologic unit called the Bearhead Rhyolite. To 
be consistent with our interpretations, in the rest of the 
paper we refer to the archeological artifacts previously 
called “Paliza Canyon obsidian” as the Bearhead 
Rhyolite obsidian. 


Geologic background of study area 


Our study area (Figs. 1, 2) encompasses both our new 
sampling localities and the Paliza Canyon sampling area 
of Nelson (1984) and Baugh and Nelson (1987). In this 
area, obsidian occurs in three types of geologic units (Fig. 
2): 1) surficial, relatively weakly consolidated tephra, 2) 
in valley-bottom alluvium, and 3) volcanic bedrock. In 
the study area, most surficial tephra is petrographically 
similar to what is now called the El Cajete Pyroclastic Beds 
of the East Fork Member of Valles Rhyolite (Gardner et 
al., 2010; Goff et al., 2011), which we simply refer to as 
El Cajete pumice, whose age is interpreted as 74.7+1.3 ka 


August 2016, Volume 38, Number 3 


53 


106°37.5' 


35°} 
45° 


106°30' 


Figure 1. Orthophoto over digital elevation model. Study area depicted by polygon. Labels show approximate location of collection localities (i.e. 
090413-2, 100915-2, etc.), geographic features, USGS quadrangle borders (in light gray), and Pueblo Revolt archaeological sites (LA numbers) 
near our study area (Liebmann 2012; Shackley 2009a, 2012b). Our precise collection localities are shown in Figure 2. Abbreviations: SR = state 
highway; FR = Santa Fe National Forest route. 


(Zimmerer et al., 2016). However, in places the El 
Cajete pumice is so thin that underlying tephras from 
older volcanic bedrock locally crop out at the surface, 
as discussed below. Geologically reworked obsidian 
is abundant in late Quaternary alluvium filling Paliza 
Canyon proper, near the boundary of the Redondo Peak 
and Bear Springs Peak Quadrangles (Fig. 2). The “Paliza 
Canyon” obsidian sampled by Baugh and Nelson in the 
1980s very likely came from late Quaternary, valley- 
bottom alluvium at this locality (Baugh and Nelson 1987; 
Tim Baugh, personal communication, 2015). 

Volcanic bedrock obsidian localities are in 6-8 Ma lavas 
and volcaniclastic sediments. These volcanic units include 
the Bearhead Rhyolite and Paliza Canyon Formation. The 
Bearhead Rhyolite consists of more than 50 felsic domes, 
flows, dikes and sills emplaced throughout the Jemez 
Mountains (Kelley et al., 2013). Exposures with obsidian 
from the Bearhead Rhyolite are located just to the north 
of the aforementioned collection areas, one of which is the 
F04-31 dome of Goff et al., (2006) (labeled in northeast 
corner of Fig. 2). The obsidian commonly comes in the 
form of marekanites (also referred to as “Apache tears”), 
which are aphyric vitreous remnants in a rhyolitic perlite 
that commonly erode into subrounded-rounded pebbles 
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(Lajédkova and Kraus 1993).” An *°Ar/*?Ar age of 7.62 
+ 0.44 Ma was obtained from an obsidian marekanite 
sample from this dome (Fig. 3; Goff et al., 2006; Kelley 
et al., 2013). Other Bearhead Rhyolite domes and flows 
in the vicinity yield ages ranging between 6.0 and 7.2 Ma 
(Justet 1996; Kempter et al., 2004). 

In contrast to the felsic Bearhead Rhyolite, the 
volcanic domes and flows in the Paliza Canyon Formation 
are all mafic to intermediate composition (Ellisor et 
al., 1996; Goff et al., 2011; Kempter et al., 2004). 
Beneath the aforementioned (often very thin) blanket 
of El Cajete pumice are irregular beds of volcaniclastic 
conglomerates and associated pumiceous deposits (Tpv in 
Fig. 2) derived primarily from andesitic to dacitic rocks 
of the Paliza Canyon Formation (Kempter et al., 2004; 
Goff and Gardner 2004; Goff et al., 2006, 2011). These 
conglomerates contain older tephras and obsidian. 


Sampling and field observations of 
obsidian sources 


Obsidian sampling efforts were undertaken in September 
of 2013 and September—October of 2015. The first trip 
was mainly a reconnaissance. In the Fall of 2015, we 
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Figure 2. Simplified geologic map of the upper Paliza Canyon area, southern Jemez Mountains, New Mexico (modified from Kempter et al., 2004, and Goff et al., 2006; labels and 
terminology from Goff et al., 2011). Qal (dark yellow) = valley bottom alluvium. Qvec (pale yellow) = El Cajete Pyroclastic Beds; moderately sorted beds of rhyolitic pyroclastic fall and 
thin pyroclastic flow deposits (74.7+1.3 ka, Zimmerer et al., 2016); locally the beds are extremely thin. Qbt (orange) = Tshirege Member, Bandelier Tuff; rhyolitic ignimbrite (1.25+0.01 
Ma, Phillips et al., 2007). Tpv (olive green) = volcaniclastic deposits, debris flows, hyper-concentrated flows, and stream deposits of the Paliza Canyon Formation. Tbh (reddish 
brown) = Bearhead Rhyolite; domes and flows of aphyric to slightly porphyritic lava (dome in northern Bear Springs Peak Quadrangle is 6.66+0.06 Ma, Kempter et al., 2004; faulted 
dome in Redondo Peak Quadrangle is 7.62+0.44 to 7.83+0.26 Ma, Goff et al., 2006). Tpa (green) = Paliza Canyon Fm, andesite, undivided; lava flows containing plagioclase and 
pyroxene (8.78+0.14 to 9.44+0.21 Ma, Justet 2003). Tpd (rose) = Paliza Canyon Fm. dacite; domes and flows of very porphyritic lava near Ruiz Peak. Tpbhd (red) = Paliza Canyon Fm. 
porphyritic biotite-hornblende dacite; dome and flows are 9.24+0.22 Ma (Justet 2003). Tpb (blue) = Paliza Canyon Fm. basalt; lava flows mostly containing visible olivine (9.45+0.07 


to 9.54+0.08 Ma, Goff et al., 2006 and Kempter et al., 2004). 


conducted sampling in order to compare the compositions 
of newly identified obsidian sources (e.g., Kempter et 
al., 2004) with the older “Paliza Canyon” source locale 
reported in Baugh and Nelson (1987). Figure 2 depicts 
our sampling localities (e.g., 090413-3); individual 
samples collected at a locality are denoted by a hyphen 
and numeral appended to the locality name (e.g., 
090413-3-2, 090413-3-3). 

On September 11, 2015, we examined geological 
localities in the northwest corner of the Bear Springs 
Peak 7.5’ USGS Quadrangle to better understand the 
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geological context of obsidian marekanites initially 
observed in mapping by Kempter et al., (2004) and 
in field reconnaissance. Because it contained obsidian 
clasts, the El Cajete pumice as mapped by Kempter et al., 
(2006) was initially considered by us as a possible source 
for the Bearhead Rhyolite obsidian. We confirmed that 
marekanites there only occurred in association with 
pyroclastic deposits, at that time thought to be El Cajete 
pumice. This pumice and accompanying obsidian was 
sampled at various localities (091115-2, 090413-2, 
090413-3, 091115-1 in Fig. 2). The largest marekanite 
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Figure 3. Photograph of Bearhead Rhyolite obsidian and dacite pumice samples by collection localities and sample numbers 
(see Figs. 1 and 2, and Table 2). Note general decrease in angularity and marekanite size from primary dome (F04-31), 
to secondary samples (locality 100915-2) in downstream alluvium in Paliza Canyon 2-3 km southwest of F04-31, locality 
091115-2 above Paliza Canyon, and locality 090413-3 farther south (see Figs. 1 and 2). Samples listed in locality 100915-2 
exclude analyzed samples in Table 2 from site 100915-2. Dacite pumice corresponds to samples 091115-2-7 and 091115-2-9. 


recovered was 55.1 mm (long axis) but most were 10-30 
mm, with a few 30-40 mm. 

The Redondo Peak 7.5’ USGS Quadrangle in Paliza 
Canyon was visited in October of 2015. There, we collected 
reworked (secondary deposit) obsidian in Paliza Canyon 
alluvium near the boundary between the Redondo Peak 
and Bear Springs Peak 7.5’ USGS Quadrangles (locality 
100915-2). This locality approximately corresponds with 
the original “Paliza Canyon” source locale (Baugh and 
Nelson 1987) and lies below and 2 to 3 km southwest 
of the unnamed dome informally called F04-31 (Figs. 
1 and 2; Goff et al., 2006; Kelley et al., 2013). In the 
Paliza Canyon alluvium, we recovered a nearly 100 
mm-long obsidian nodule and found abundant (100s 
per Sm’) obsidian marekanites (see sample 100915-2 in 
Figure 3). Many of the marekanites were around 50 mm 
in diameter (Fig. 3). The purpose of the pumice samples 
was to determine if their composition was rhyolitic (and 
thus Bearhead Rhyolite) or dacitic (and thus part of the 
Paliza Canyon Formation). 


XRF analyses 


We used XRF to determine the chemical composition of 
these obsidian samples and tephra from the Paliza Canyon 
Formation. XRF analytical instrument methodology is 
detailed in Appendix 1 (see also Shackley 2005). Major 
oxide results are presented in Table 1 and Figure 4, and 
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Table 2 lists selected trace element concentrations. In 
Table 3, we compare the mean and central tendency 
for the trace element data of our samples with those of 
Baugh and Nelson’s (1987) original study. 

Of primary importance to our study is the 
compositional similarity of the F04-31 rhyolite dome 
with alluvial obsidian samples collected at the “Paliza 
Canyon” source locality. This similarity is especially 
evident in the overlap of trace element concentrations 
(Table 2) between samples 100915-2-n (“Paliza 
Canyon” source locality) and F04—31-n (Bearhead 
Rhyolite dome). The mid-Z incompatible elements 
(Rb, Sr, Zr) most frequently used in geoarchaeological 
studies are generally similar, with minor differences in 
mean elemental concentrations (Table 3) attributable to 
changes in analytical technology since the 1980s (see 
Shackley 2011). Specifically, modern Si(Li) and SDD 
EDXRF detectors are more sensitive than the earlier 
WXRF instruments used for Baugh and Nelson’s 
(1987) four samples (Shackley 2011). Figures 5 and 6 
nicely illustrate this overlap between the two localities 
(100915-2 and F04-31) for Nb, Y, Ba concentrations. 
Thus, we interpret that the FO4-31 dome is the source 
for the downstream alluvial obsidian clasts found in the 
“Paliza Canyon” source locality of Baugh and Nelson 
(1987). Note the relatively high Sr and Ba values of 
samples correlative to the Bearhead Rhyolite (brown, 
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TABLE 1. Major and minor oxide data for pumice, obsidian, and lava rock samples. 


Sample SiO, Al,O, CaO Fe,0, K,O MgO MnO Na,O TiO, = 
091115-2-7! 67.96 17.06 1,29 5.74 3.81 1.07 0.13 1:91 0.78 99.74 
091115-2-8! 64.06 19.70 1.49 7.002 3.82 1.16 0.08 1.66 0.71 99.69 
091115-2-9! 68.61 15.73 1.65 5:93 3.72 1.35 0.08 1.83 0.72 99.62 
090413-3-22 74.67 13.10 0.89 1.48 5.42 0.21 0.13 3.37 0.28 99.55 
Banco Bonito F82-7? 73.9 13.6 159 1.86 4.15 0.64 0.05 3.87 0.29 99.95 
EI Cajete F82-64 72.8 14.0 1.87 2.11 4.16 0.78 0.05 3.84 0.33 99.94 
Paliza Cyn dacite° 68.99 13.35 2.62 3.84 2.96 0.77 0.09 3.97 0.74 99.75 
RGM-1 (this study) 75.68 12.48 1.30 1.81 4.55 <0.1 0.04 Salt 0.20 99.82 
RGM-1 USGS 73.4 13:7 1.15 1.86 4.30 0.28 0.036 4.07 0.27 99.06 
recommended 


‘Pumice samples associated with obsidian localities formerly mapped 
(Kempter et al., 2004). 

? Obsidian sample from an area formerly mapped as El Cajete pumice 
3from Gardner et al. (1986) 

“from F. Goff (unpub. data) 

5 from Ellisor et al. (1996) 

Data not normalized to USGS recommended values. 


as El Cajete pumice above and east of Paliza Canyon 


above and east of Paliza Canyon (Kempter et al., 2004) 
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Figure 4. TAS plot (Le Bas et al., 1986) of selected samples in this study from Table 1. The elemental and major oxide data in Table 1 and 


this TAS plot demonstrate that some of the pumice in our study area is 


dacitic and probably from the volcaniclastic member of the Paliza 


Canyon Formation (Neogene), rather than rhyolitic pumice of the Quaternary El Cajete pumice (Qvec) (c.f. Kempter et al., 2004; see Figure 
2 herein). While a seemingly minor issue, it does point to the mixing of obsidian sources in the southern Jemez Mountains and presents 


another complication for geoarchaeological investigations. 


periwinkle gray, light green, purple, tan text in Table 
2) compared to other archaeologically relevant rhyolites 
in the Jemez Mountains: Valles, Cerro Toledo, El 
Rechuelos, and Canovas Canyon (lower part of Table 2). 

Obsidians collected from areas mapped as El Cajete 
pumice in the Bear Springs Peak Quadrangle, southeast 
of Paliza Canyon, have a composition compatible with 
the Bearhead Rhyolite. Comparison of trace element 
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concentrations of these samples indicate overlap with 
the obsidian from the F04-31 Bearhead Rhyolite dome 
(Figures 5 and 6) but notable dissimilarity with El 
Cajete pumice and Banco Bonito, the latter two being 
notably similar (Table 2). Samples of pumice from these 
obsidian localities also generally do not chemically 
match the El Cajete pumice. One obsidian sample had 
a similar chemistry to the Canovas Canyon Rhyolite 
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Figure 5. Nb, Y, Ba three-dimensional plot of Bearhead Rhyolite obsidian data in Table 
2 and the one Canovas Canyon Rhyolite obsidian marekanite from collection locality 
091115-1 (see Fig. 2). Color-coding of collection localities follows that of Table 2. 
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Figure 6. Nb versus Y bivariate plot of Bearhead Rhyolite obsidian data from Table 
2. Color-coding of collection localities is from Table 2. Note that the primary source 
obsidian from the F04—-31 Bearhead Rhyolite dome (light blue squares) occupies 
the same space as secondary deposit samples from Paliza Canyon (100915-2) and 
collection localities south of Paliza Canyon (see Fig. 2). 
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TABLE 2. Selected trace element concentrations for obsidian and pumice samples 


Sample* Zn Rb Sr Y Zr Nb Ba Pb Th 


090413-3-1 56 109 of we 128 34 1741 24 1 
090413-3-11 61 100 87 DS il 36 1612 ty 17 
090413-3-12 58 101 86 22, 120 34 1665 19 ip 
090413-3-14 32) 103 86 28 128 35) 1685 19 16 
090413-3-18 78 105 87 24 122 32 1672 Dil 12 
090413-3-19 101 Oy 85 23 itil) 30 1689 18 13 
090413-3-2 a7 oS) 85 DS 124 38 1607 15 iis 
090413-3-20 81 102 82 24 120 31 1577 22 10 
090413-3-21 87 95 81 23 lS 27, 1707 17. ily 
090413-3-23 89 100 85 DS 119 27 1656 UP 16 
090413-3-24 al 93 79 22 128 33 1606 18 10 
090413-3-3 47 99 86 26 120 32 1640 23 iil 
090413-3-4 46 OY) 86 22 136 88) 1445 18 12 
090413-3-5 56 105 89 20 138 ay 1706 2p) iQ) 
090413-3-6 Auy 101 86 20 1S 34 1680 19 i) 
090413-3-7 53 105 93 oA 128 YY) 1560 aD 7 

090413-3-8 55) 102 89 24 sy 29 1609 20 13 
090413-3-9 ae 96 85 24 130 31 1661 22) ily 
091115-1-1 EPS) 93: 85 26 140 aril 1651 18 17 
091115-1-3 43 100 86 26 130 ill 1997 19 8 

091115-1-5 56 111 Oil 26 135 34 1630 DS 18 
091115-1-6 49 102 86 NY 129 35 1699 19 20 
091115-1-7 59 98 86 NY 132 33 1673 19 14 
091115-1-8 ott 101 88 25 129 of 1668 20 13 
091115-1-9 ot 105 89 28 137 36 1601 20 US 


100915-2-1 (“Paliza Canyon” source area) 49 104 90 AY, 135 36 1661 21 14 
100915-2-10 (“Paliza Canyon” source area) 50 95) 84 25) 7; 55) 1879 18 10 
100915-2-2 (“Paliza Canyon” source area) 54 108 94 Diy i333} 36 1693 19 11 
100915-2-3 (“Paliza Canyon” source area) 42 94 80 28 121 Sil 1642 16 11 
100915-2-4 (“Paliza Canyon” source area) Dill 105 OD) 27 140 ayy 1778 20 14 
100915-2-5 (“Paliza Canyon” source area) 45 98 88 24 2D) 32 1639 20 19) 
100915-2-6 (“Paliza Canyon” source area) 44 108 3} 22 33 32 193; AA 13 
100915-2-7 (“Paliza Canyon” source area) 44 OY) OE) 23} 128 34 1643 22) 14 
100915-2-8 (“Paliza Canyon” source area) 41 90 T 23} 22} 30 1646 16 11 
100915-2-9 (“Paliza Canyon” source area) 50 102 94 27 32! Sil 1686 Dl 15 


Table 2 continued next page. 
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Table 2 continued from previous page. 


Sample* Zn Rb Sr Y Zr Nb Ba Pb Th 
F04-31-1 (Bearhead Rhyolite Dome) 45 91 79 AS) 123 32 1883 15 18 
F04-31-2 (Bearhead Rhyolite Dome) 46 104 90 2S) 142 B2 1766 D?) 16 
F04-31-3 (Bearhead Rhyolite Dome) 44 100 85 aT 132 28) 1655 Ail ily 
Valles Rhyolite (Cerro del Medio)! 26° 160 10 43 172 54 35 27.7% 17.5° 
Cerro Toledo Rhyolite! Le. 207 5 63 183 98 23 mt 23} 
EI Rechuelos Rhyolite! e. 152 9 23 77 47 24 me 16.75 
Canovas Canyon Rhyolite! 217 116 43 21 108 53 352 mt. 129 
El Cajete pumice® 32? 136 185 27 125 36 616 19 18 
Banco Bonito Rhyolite obsidian® 30% 143 192 28 136 39 596 20 18.6 
RGM-1 (this study) 39 151 104 27 218 8 806 24 18 
RGM-1 USGS recommended 32 150 110 25 220 8.9 810 24 15 


*39 Bearhead Rhyolite obsidian samples were collected from areas previously mapped as El Cajete pumice (brown, periwinkle gray, light green, and purple text for 090413-2-n, 
090413-3-n, 091115-1-n, and 091115-2-n respectively), 10 samples of alluvial obsidian from the “Paliza Canyon” source area (tan text), and 3 samples of obsidian from the F04-31 
Bearhead Rhyolite dome (light blue text; Goff et al., 2006). Comparative samples include: El Cajete pumice, one sample each from known artifact quality obsidian sources from the 
Jemez Mountains, and the USGS RGM-1 rhyolite standard (see Shackley 2005). Note that color-coded localities correspond with colors plotted in elemental plots (Figures 5 and 6). 
‘Mean concentrations mainly from unit Qvdmw of Gardner et al. (2007); see Appendix in Shackley 2005. 


°Zn concentration from F. Goff, unpub. data. 

3Zn concentration from Gardner et al. (1986). 

“Th concentration from Stix et al. (1988). 

5Th concentration from Loeffler et al. (1988). 

6Zn, Pb and Th concentrations from Gardner et al. (2007); obsidian unit Qvdmw. 
Th concentration from Gardner et al. (1986). 


Except for Zn, the El Cajete pumice and Banco Bonito obsidian values are from Self et al. (1988). 


n.r. = not reported. 


TABLE 3. Mean and central tendency for Bearhead Rhyolite obsidian from data in Table 2 and 
mean data from the four “Paliza Canyon” samples originally reported by Baugh and Nelson (1987). 


Element N Minimum Maximum Mean‘ oO 

Zn 61 41 145 60 (n.r.)? 22:5 
Rb 61 87 111 100 (100.8) 5:3 
Sr 61 76 99 87 (86.2) dy 
Y 61 20 29 25 (8.5) 2.2 
Zr 61 112 142 128 (124) 7.1 
Nb 61 27 41 33 (18.1) 2.8 
Ba 61 1445 1883 1660 (1352) 7a 
Pb 61 15 25 20 (n.r.) 2.1 
Th 61 7 21 14 (n.r.) 3.1 


‘ Baugh and Nelson (1987) mean values reported in parentheses. 
? nr. = not reported by Baugh and Nelson (1987). 


(sample 091115-1-10; Table 2, Fig. 5). The remaining 
tephra samples were dacitic (091115-2-7, 091115-2-8, 
and 091115-2-9 samples in Table 1 and Fig. 4) and likely 
part of the Paliza Canyon Formation, as discussed below. 


Discussion 


Sources for the formally known “Paliza 
Canyon” obsidian 


The geochemistry in the preceding section demonstrates 
that the “Paliza Canyon obsidian” in the archaeological 
vernacular compares very well with Bearhead Rhyolite 
obsidian, consistent with geologic relations near the 
obsidian source locality of Nelson (1984) and Baugh and 
Nelson (1987). We interpret that the primary sources of 
this obsidian are domes of Bearhead Rhyolite near the 
head of Paliza Canyon in the southern Jemez Mountains, 
one of which is the FO4—31 dome of Goff et al., (2006). 
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Small but fresh obsidian fragments (marekanites or 
Apache tears) are also associated with the Peralta Tuff 
Member of the Bearhead Rhyolite (ca. 7 Ma) in the 
southeastern Jemez Mountains (F. Goff, unpublished 
data). Whether the Peralta Tuff material is artifact quality 
obsidian from an unknown locale is not yet clear. Other 
geologic source areas near Paliza Canyon are discredited 
for reasons enumerated below. 

Geologists working in the Jemez Mountains have 
noted that obsidian is abundant in the upper part of the 
El Cajete pumice that outcrops near all the collection 
localities (references above, especially Goff et al., 2006; 
Kempter et al., 2004). In fact, we initially hypothesized 
that the El] Cajete pumice was the source of what we now 
call Bearhead Rhyolite obsidian. However, the obsidians 
from the sampled tephras southeast of Paliza Canyon 
in the Bear Springs Peak Quadrangle chemically match 
the Bearhead Rhyolite, as noted above. Furthermore, 
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the associated tephra are dacitic in composition, based 
on hand lens examination plus major and minor oxide 
chemistry data (Table 1; TAS plot in Fig. 4). The dacitic 
tephras are probably from dome eruptions associated 
with the Paliza Canyon Formation, such as unit Tpbhd 
to the north (Fig. 2), or pumice from various Paliza 
Canyon Formation domes reworked into volcaniclastic 
gravels (Tpv) of the Keres Group (see Goff et al., 
2006; Kempter et al., 2004). These older Keres Group 
tephras and gravels appear to be overlain by a thin and 
discontinuous cover of El Cajete pumice, at least in our 
sampled area above and southeast of Paliza Canyon. 
For example, tephra and reworked obsidian from unit 
Tpv underlies a very thin blanket of El Cajete pumice at 
sampling locality 091115-1. 

Another possible source is the Canovas Canyon 
Formation, which locally has obsidian. However, 
Canovas Canyon Formation obsidian is readily 
differentiated from Bearhead Rhyolite obsidian by 
their lower Sr and Ba concentrations and, to a lesser 
extent, differences in Y and Nb composition (Table 
2; Fig. 5). One of marekanite samples from collection 
locality 090413-3/091115-1 (same locality, two 
separate collection dates), in what we thought was El 
Cajete pumice, had similar chemistry to the Canovas 
Canyon Formation (Table 2, sample 091115-1- 
10)-although all the other collected obsidian at this 
locality was Bearhead Rhyolite obsidian (Fig. 5). Like 
the aforementioned Paliza Canyon dacitic tephras, 
the marekanite represented by sample 091115-1-10 
is interpreted to have been reworked into the Keres 
Group from a Canovas Canyon Formation source. The 
Canovas Canyon Rhyolite is mapped just southeast of 
our study area, so this is not unexpected (see Fig. 2 and 
Kempter et al., 2004). 

Despite their aforementioned geochemical differences, 
the Bearhead Rhyolite obsidian and Canovas Canyon 
Rhyolite obsidian are characterized by the highest Sr and 
Ba and lowest Rb of all the Jemez Mountain sources (Table 
2), indeed all the archaeological obsidian from Jemez 
Lineament obsidian sources (Mount Taylor and Jemez 
Mountains in Shackley 2005). Both Bearhead Rhyolite 
obsidian (ca. 7 Ma) and Canovas Canyon Rhyolite (Bear 
Springs Peak) obsidian (ca. 8.85-9.7 Ma) have been eroding 
down into the Rio Grande for millions of years (Goff et al., 
2006; Kempter et al. 2004). While representing relatively 
small percentages of obsidian provenance assemblages in 
North American Southwestern archaeological sites, they 
nevertheless are geochemically different than the more 
common Jemez Mountains obsidian sources such as Cerro 
Toledo Rhyolite, El Rechuelos Rhyolite, and Valles Rhyolite 
(Cerro del Medio). 


Archeological Implications 


The largest Bearhead Rhyolite marekanite recovered was 100 
mmin largest dimension but most were 10-30 mm, witha few 
30-40 mm. Given that recycled Bearhead Rhyolite obsidian 
recovered at Tijeras Wash (south of Albuquerque, more than 
150 stream km south of the southern Jemez Mountains) are 
as much as 3 mm in largest dimension, obsidian from the 
primary sources must have been significantly larger than 
the marekanites recovered in this study. In October 2015, 
we recovered a marekanite nearly 100 mm in diameter 
(sample 100915-2-1 in Table 2 and Fig. 3) in alluvium 
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filling Paliza Canyon just downstream from the F04-31 
Bearhead Rhyolite dome, certainly large enough to supply 
samples downstream through fluvial transport along the 
Rio Grande. As a media for tool production, this smoky 
gray to nearly transparent obsidian is equal in quality to 
other Jemez Mountains sources, as the Ancestral Puebloans 
found during the Pueblo Revolt period. 

Pueblo Revolt period knappers in the area used this 
obsidian due to proximity, but it was also used as a raw 
material throughout prehistory in the region (Shackley 
2014a, 2014b, 2015; see below). Other than archaeological 
contexts in the immediate Jemez Mountains region, one 
could argue that this Bearhead Rhyolite obsidian is a minor 
component of archaeological obsidian raw material used in 
the Southwest. These “minor” sources do have potential to 
illuminate past human behavior and can be as important 
as major sources, as has been argued elsewhere (Shackley 
2009b, 2009c). A recent example from northern New 
Mexico proves illuminating. 

At the Late Archaic Early Agricultural site of San Luis 
de Cabezén northwest of Cuba, New Mexico (LA 110946), 
the obsidian assemblage was dominated by artifacts, 
mainly debitage (waste flakes) produced from Cerro Toledo 
Rhyolite (64.4%) and Valles Rhyolite (Cerro del Medio; 
25%), the sources of which are about a day’s walk away in 
the Jemez Mountains to the east (n=104 obsidian artifacts; 
Shackley 2014b). What is notable is that many of the 
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El Rechuelos 
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(Cerro del Medio) Bear Springs Peak 
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Figure 7. Late Archaic obsidian projectile points and source provenance 
from San Luis de Cabezén (LA 110946) produced from Bearhead 
Rhyolite (previously called Paliza Canyon) and Canovas Canyon Rhyolite 
(Bear Springs Peak) obsidian, and one produced from Valles Rhyolite 
and El Rechuelos Rhyolite obsidian (from Shackley 2014b). Projectile 
point 5749 is nearly 35 mm in length produced from a Bearhead Rhyolite 
obsidian core at least 40 mm in diameter. Note: spherulite in sample 
5745 typical of Cerro del Medio obsidian and generally not present in 
the Keres Group obsidians. 
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obsidian projectile points recovered (Fig. 7) were produced 
from Keres Group sources that are in the southern portion 
of the Jemez Mountains and nearer to this site. These 
Keres Group sources include Bearhead Rhyolite (previously 
called “Paliza Canyon”; 4.8% of total obsidian), Canovas 
Canyon Rhyolite (3.8% of total obsidian), the only artifact 
produced from Polvadera Group El Rechuelos on the north 
side of Valles Caldera, and one projectile point produced 
from Tewa Group Valles Rhyolite, although Valles Rhyolite 
debitage comprised 25% of the overall assemblage (Fig. 
7; Shackley 2014b). The Keres Group sources are small 
Neogene marekanite sources, but provide excellent media 
for flaked stone tool production and are located nearer 
to the Cabezén site than most of the other Jemez sources. 
Additionally, projectile points produced from Bearhead 
Rhyolite obsidian at San Luis de Cabez6n were near 50 mm 
in length, and so had to be produced from cores larger than 
50 mm in size—within the range of marekanites recovered in 
Paliza Canyon just below the F04-31 dome. 

The above discussion suggests that the raw material 
at San Luis de Cabez6n was procured from the primary 
dome or nearby in the southern Jemez Mountains rather 
than from secondary deposits to the south. Therefore, the 
presence of “minor” obsidian sources in this case elucidates 
prehistoric hunting and raw material procurement behavior 
that might not be obvious looking at the major sources in the 
Jemez Mountains alone. Indeed, understanding the primary 
location of a raw material source is crucial 1 in more precisely 
defining behavior, even if the source is a “minor” source. 

The unity of geology and archaeology, a long-standing 
relationship, can offer clarity to both disciplines. An 
entire Geological Society of America volume was devoted 
to this relationship (Stein and Linse 1993). The Bearhead 


Rhyolite obsidian project is a prime example of this 
“symbiotic” relationship, and in part shows the differences 
in scientific scale between these closely related disciplines. 
Our discovery of dacitic tephra (reworked from the Paliza 
Canyon Formation) in what has previously been mapped 
as El Cajete Rhyolite is mainly due to a finer-grained 
approach common in archaeology. Mapping with this 
finer-grained approach, particularly in archaeological sites, 
is often required to be at the millimeter level, a scale rare in 
geological mapping. This archaeological perspective is often 
brought into geoarchaeological investigations and, as seen 
here, has ramifications for geological interpretation. 
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Appendix 


The ThermoFisher/Scientific QUANT’X 
EDXRF Spectrometer methods at the 
Geological XRF Laboratory 


This appendix presents analytical methodology for the 
EDXRF Spectrometer used in this study (modified slight- 
ly from http://swxrflab.net/anlysis.htm). All archaeolog- 
ical samples are analyzed whole. The results presented here 
are quantitative in that they are derived from “filtered” inten- 
sity values ratioed to the appropriate x-ray continuum regions 
through a least squares fitting formula rather than plotting the 
proportions of the net intensities in a ternary system (McCar- 
thy and Schamber 1981; Schamber 1977). Or more essentially, 
these data, through the analysis of international rock stan- 
dards, allow for inter-instrument comparison with a predict- 
able degree of certainty (Hampel 1984). 

The spectrometer is equipped with: (1) a peltier air cooled 
solid-state Si(Li) X-ray detector having an ultra ultra-high- 
flux end window bremsstrahlung rhodium (Rh) x-ray target 
with a 125 micron beryllium (Be) window and an 8.8 mm 
collimator (3.5 mm for samples < 10mm diameter), and (2) 
x-ray generator that operates from 4-50 kV/0.02-1.0 mA 
at 0.02 increments using an IBM PC based microprocessor 
and WinTracelM 7.0 software. The spectrometer is equipped 
with a 2001 min-1 Edwards vacuum pump for the analysis of 
elements below titanium (Ti). Data is acquired with a pulse 
processor and analog to digital converter. This is a significant 
improvement in analytical speed and efficiency beyond the 
former Spectrace 5000 and QuanX analog systems (see Davis 
et al., 2011; Shackley 2005). 


Trace Element Analyses 


For Ti-Nb, Pb, Th elements the mid-Z(b) group condition is 
used operating the x-ray tube at 30 kV. The x-ray tube utilizes a 
0.05 mm (medium) Pd primary beam filter in an air path at 200 


seconds livetime to generate x-ray intensity Kal-line data for 
elements titanium (Ti), manganese (Mn), iron (as FeT), cobalt 
(Co), nickel (Ni), copper (Cu), zinc (Zn), gallium (Ga), rubidi- 
um (Rb), strontium (Sr), yttrium (Y), zirconium (Zr), niobium 
(Nb), lead (Pb), and thorium (Th). Not all these elements are 
reported since their values in many volcanic rocks are very low 
and below detection limits. Trace element intensities were con- 
verted to concentration estimates by employing a least-squares 
or quadratic calibration line ratioed to the Compton scatter 
established for each element from the analysis of international 
rock standards certified by the National Institute of Standards 
and Technology (NIST), the US. Geological Survey (USGS), 
Canadian Centre for Mineral and Energy Technology, and 
the Centre de Recherches Pétrographiques et Géochimiques in 
France (Govindaraju 1994). Line fitting is linear (XML) for all 
elements. When barium (Ba) data are acquired, the Rh tube is 
operated at 50 kV and .25-1.0 mA (automated) in an air path 
at 200 seconds livetime to generate x-ray intensity Ka1-line 
data, through a 0.630 mm Cu (thick) filter ratioed to a portion 
of the bremsstrahlung region (see Davis et al., 2011). Further 
details concerning the petrological choice of these elements 
in Southwest obsidians is available in Shackley (2005; also 
Mahood and Stimac 1991; and Hughes and Smith 1993). A 
suite of 16 specific U.S. Geological Survey standards are used 
for the best fit regression calibration for elements Ti- Nb, Ba, 
Pb, and Th: G-2 (basalt), AGV-2 (andesite), GSP-2 (granodi- 
orite), SY-2 (syenite), BHVO-2 (hawaiite), STM-1 (syenite), 
QLO-1 (quartz latite), RGM-1 (obsidian), W-2 (diabase), 
BIR-1 (basalt), SDC-1 (mica schist), BCR-2 (basalt), TLM-1 
(tonalite), SCO-1 (shale), NOD-A-1 (manganese), NOD-P-1 
(manganese). Four other standards include: National Institute 
of Standards and Technology NIST-278 (obsidian), BR-E 
(basalt) from the Centre de Recherches Pétrographiques et 
Géochimiques in France, and JR-1 and JR-2 (obsidian) from 
the Geological Survey of Japan (Govindaraju 1994). 

The data from the WinTrace software are translated 
directly into Excel for Windows software for manipulation and 
into SPSS for Windows for statistical analyses when necessary. 


TABLE 1. X-ray fluorescence concentrations for selected trace elements from obsidian standard RGM-1* 


SAMPLE Ti Mn Fe Rb Sr Y Zr Nb Ba Pb Th 

RGM-1 

(Govindaraju 1600! 279 12998 149 108 A) 219 8.9 807 24 15.1 

1994) 

RGM-1 (USGS 1618+120 279+50 12998+210 150+8 110+10 257 220+20 8.940.6  810+46 24+3 15+1.3 

recommended)! 

GeoRem values*  1510- 282 12788- 142-164 96.73-116 21.6-25.1 173-258 8.37-13 791-881 18-284 1416.3 
1940 14263 

RGM-1, pressed 1523249 294413 13723430 149+2 108+2 2542 219+2 9+2.0 8043 21+1.6 17+3.2 

powder (this 

study, n=99) 

RGM-1, flake 1568+44 3114+11 13306+33 153+2 113+2 2541.5 230+4 9+2 942414 23+1.5 1543.4 

from original 

USGS boulder 

(n=12) 


* Pressed powder and whole rock elemental analysis of USGS RGM-1 obsidian standard measured on the ThermoScientific Quant'X, which are compared to USGS, 
Govindaraju (1994), and GeoRem recommended values. USGS RGM-1 pressed powder pellet (n=99 runs) and whole rock flake are from original USGS boulder (n=12). 
+ values represent first standard deviation computations for the group of measurements. All values are in parts per million (ppm) as reported in Govindaraju (1994), 
USGS, and this study. RGM-1 is a U.S. Geological Survey obsidian standard obtained from Glass Mountain, Medicine Lake Highlands Volcanic Field, northern California. 


‘Ti, Mn, Fe calculated to ppm from wt. percent from USGS data using MURR’s oxide to element and element to oxide multiplier table (from Glascock 1991): http://swxrflab. 


net/MURRoxidetoelementtable.jpg 
? information value 
3n=17 
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In order to evaluate these quantitative determinations, machine 
data were compared to measurements of known standards 
during each run (Table 1). RGM-1 is analyzed during each 
sample run (n=1-19) for obsidian artifacts to check machine 
calibration. Other appropriate standards from the above list 
are used for other volcanic rocks. 


Detection limits 


Detection limits for EDXRF have been calculated at 6 = 
at this web site: http://swxrflab.net/detectionlimits.htm. 
Recently, a number of scholars have questioned the validity 
of using pressed powder pellets of international standards for 
empirical calibration and data checking. This concern, called 
a “matrix issue,” questions the potential differing analytical 
results between pressed powder and whole rock samples (Jeff 
Ferguson, personal communication, 2010; cf. Shackley 2011). 
The potential problem was tested in detail in Kathleen Davis 
et. al’s, (1998, reprinted 2011) study of EDXRF on obsidian 
and found to not be an issue (see also Shackley and Hampel 
1992). Hermes and Ritchie (1997) derived similar conclusions 
using EDXRF with archaeological felsites. Table 1 exhibits the 
analysis of one of my pressed powder pellets of USGS RGM-1 
and a flake from the same USGS collected boulder sent by 
Steve Wilson of the USGS. As one can see, the variability is 
within 1% or less in most cases, and some of this variability 
is inherent in the native variability within the single 200 kg 
boulder collected at Glass Mountain by USGS, not because of 
differing matrices. This issue is well examined by Ron Jenkins 
and seems to be misunderstood in the discipline (c.f. Jenkins 
1999:167-172). While it is true that his discussion hinges on 
the assumption of homogeneity, one must assume that for 
purposes of non-destructive EDXRF that pressed powder 
standards are heterogeneous while obsidian is homogeneous. 
This is, indeed, not the case for obsidian studies in general. 
Jenkins explains: 


“Primary absorption occurs because all atoms of the 
specimen matrix will absorb photons from the pri- 
mary source. Since there is a competition for these 
primary photons by the atoms making up the spec- 
imen [pressed powder versus glass], the intensity/ 
wavelength distribution of these photons available 
for the excitation of a given analyte element may be 
modified by other matrix effects (Jenkins 1999:168; 
emphasis mine).” 


In the case of pressed powder pellets, the “competition” 
involves air molecules locked in the powder matrix, which 
is not analyzed in obsidian studies and of course is not 
present in volcanic glass, particularly with XRF and ICP- 
MS technology. While it is true that the borate mix used in 
the pressed powder can yield some boron, this element is 
undetectable by XRF and so light (Z=5) as to not cause peak 
overlap, especially when using any tube filtering. Perhaps 
more important is the issue of infinite thickness, particularly 
when Z>51 and high tube voltages are required. In this 
case, some of the x-rays are radiated through the sample, 
and ratioing to the Compton or bremsstrahlung regions can 
yield complex results, as can be seen in Table 1 for Ba and 
evidently a problem for other analysts as well, as seen by the 
GeoRem values in Table 1. In this case, careful calibration 
can mitigate some of this error. Parenthetically, only recently 
have PXRF products been able to reach 50 kV in order to 
adequately move electrons out of orbit for those atoms Z>51, 
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despite opinion to the contrary (see Jenkins 1999:9-12; 
Speakman and Shackley 2013). 

As I’ve noted in print and in discussions with many of 
my colleagues, not using international standards, as is often 
seen in PXRF studies, severely hinders the potential for 
establishing validity in obsidian research. Unless international 
standards are used, it is nearly impossible to compare results 
across laboratories (Shackley 2010, 2011; Speakman and 
Shackley 2013). Therefore, using the “matrix” issue as an 
excuse is unacceptable. 


Major Oxide Analyses 


Analysis of the major oxides of Si, Al, Ca, Fe, K, Mg, Mn, Na, 
and Ti is performed under the multiple conditions elucidated 
in Table 2. This fundamental parameter analysis (theoretical 
with standards), while not as accurate as destructive analyses 
(pressed powder and fusion disks), is usually within a few 
percent of actual based on the analysis of USGS RGM-1 
obsidian standard (see also Shackley 2011). The fundamental 
parameters (theoretical) method is run under conditions 
commensurate with the elements of interest and calibrated 
with 11 USGS standards (RGM-1, rhyolite; AGV-2, andesite; 
BHVO-1, hawaiite; BIR-1, basalt; G-2, granite; GSP-2, 
granodiorite; BCR-2, basalt; W-2, diabase; QLO-1, quartz 
latite; STM-1, syenite), and one Japanese Geological Survey 
rhyolite standard (JR-1). See Lundblad et al., (2011) for an 
alternative set of conditions and methods for oxide analyses. 


TABLE 2. Conditions of Fundamental Parameter Analysis! 


Low Za (Na, Mg, Al, Si, P) 


Voltage 6 kV Current Auto? 
Livetime 100 seconds Counts Limit 0 

Filter No Filter Atmosphere Vacuum 
Maximum 10 keV Count Rate Low 
Energy 

Mid Zb (K, Ca, Ti, V, Cr, Mn, Fe) 

Voltage 32 kV Current Auto 
Livetime 100 seconds Counts Limit 0 

Filter Pd (0.06 mm) Atmosphere Vacuum 
Maximum 40 keV Count Rate Medium 
Energy 

High Zb (Sn, Sb, Ba, Ag, Cd) 

Voltage 50 kV Current Auto 
Livetime 100 seconds Counts Limit 0 

Filter Cu (0.559 mm) Atmosphere Vacuum 
Maximum 40 keV Count Rate High 
Energy 

Low Zb (S, Cl, K, Ca) 

Voltage 8 kV Current Auto 
Livetime 100 seconds Counts Limit 0 

Filter Cellulose (0.06 mm) Atmosphere Vacuum 
Maximum 10 keV Count Rate Low 
Energy 


‘Multiple conditions designed to ameliorate peak overlap identified with digital 
filter background removal, least squares empirical peak deconvolution, gross 
peak intensities and net peak intensities above background. 

? Current is set automatically based on the mass absorption coefficient. 
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